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// a sequence of digitally on-off modulated optical pulses is injected 
into a dielectric waveguide, these pulses may begin to overlap after a 
sufficient distance of propagation because of material dispersion and/ or 
group delay spreading. In general, the pulses will not add linearly in 
power, which can complicate the problem of equalization of the square-law 
(power) detected overlapping output pulses at baseband. This paper 
illustrates important situations in which the guide may be treated as 
"pseudo-linear" in power, meaning that the detected guide output pulses 
appear to add linearly. 

I. INTRODUCTION 

If a single pulse of optical energy propagates along a dielectric 
waveguide, pulse broadening can occur for one or more of the following 
reasons : material dispersion, individual mode waveguide dispersion, or 
differences in the group delays of different guide modes. In addition, 
the pulse shape may become only statistically defined because of 
random mode coupling and/or statistical fluctuations of the optical 
source. 

If a sequence of digitally on-off modulated pulses is injected into a 
dielectric waveguide, those pulses may begin to overlap after a sufficient 
distance of propagation. In general, the optical powers in the pulses 
will not add in a linear manner.* On the other hand, as will be shown 
below, the guide may be pseudo-linear in power. That is, for the 
purpose of processing the power received at the output end of the 



t The fiber is a medium which is linear in E field propagation (from Maxwell's 
equations). It is usually excited by a power-modulated source, and its output field is 
detected by a square-law (power sensitive) device. Even if the input pulses are 
separate, the response of the square-law detector will in general contain cross terms 
resulting from the overlap in the output pulses. 
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guide — in order to make decisions as to whether or not each pulse is on 
or off — we may be able to treat the individual overlapping output 
pulses as if they added linearly. 

If the output power pulses could be considered to add linearly, then, 
after detection, the resulting current pulses could be separated with a 
linear equalizer provided the shapes of the pulses are well denned and 
identical from pulse to pulse. In this paper we consider a number of 
interesting cases in which the guide can be treated as if it were linear 
in power (output pulses add linearly) and in which the pulses at the 
output assume a well-defined shape in spite of random mode coupling 
and/or source fluctuations. 

II. CASE i: MULTIMODE GUIDE, MODE-LOCKED SOURCE, NO MODE 
COUPLING 

The easiest case to visualize in which the guide appears linear in 
power with overlapping output pulses is that of a multimode guide 
propagating pulses derived from a mode-locked laser operating in a 
single spatial mode. It is assumed that there is no mode conversion 
and that group delay differences among the modes dominate pulse 
spreading. 

Even though the laser puts out a well-defined spatial mode, it may 
be very difficult to match this to a given fiber mode so that only one 
fiber mode is excited. We assume that a number of fiber modes are 
excited by the pulses from the mode-locked laser. The assumption of a 
mode-locked laser implies that the optical bandwidth being used is 
small so that material and waveguide dispersion can be neglected, and, 
in addition, no random fluctuations are present from beating of un- 
locked source modes. The sequence of nonoverlapping pulses from the 
laser exciting the guide input will produce identical sequences of non- 
overlapping pulses in each mode at the guide output. However, the 
sequences at the output in the various modes will have relative time 
delays because of the differences in the group delay per unit length 
associated with the various modes (see Fig. 1). When the fiber output 
falls upon a detector, the current produced (neglecting shot noise) is 
proportional to the sum of the powers in all the modes. The sum of the 
powers in all the modes resulting from a single input pulse is shown in 
Fig. 1. Since the output pulses in a single mode do not overlap and 
since the detector linearly adds the powers of the various modes, the 
total detected current will be a sum of pulses modulated on and off, 
each of which looks like the response to a single input pulse. Thus, the 
output power produces a detected current which is a filtered version of 
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Fig. 1 — Multimode propagation without coupling. 



the input power. That is, the sequence of nonoverlapping input pulses 
produces a detected current at the output, which is a sequence of over- 
lapping pulses that add linearly. For the purposes of processing this 
current, the guide can therefore be considered linear in power or linear 
at baseband. 



III. CASE II : INCOHERENT SOURCE, MULTIMODE GUIDE, NO MODE 
COUPLING 

In this example, we assume that a pulse modulated incoherent 
source excites one or more modes of a multimode waveguide with no 
mode coupling. We show that the received output power in each mode 
can be treated as a linearly filtered version of the input power, i.e., that 
the sequence of nonoverlapping input pulses produces a sequence of 
output pulses in each mode which add up linearly, whether they overlap 
(because of material dispersion) or not. Since the detector produces a 
current which is proportional to the sum of the powers in each mode, 
the total current will also consist of a sequence of pulses which add 
linearly, whether they overlap or not. We show that, in order for this 
effective linearity in power to be valid, it is necessary that the source 
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bandwidth be sufficiently large compared to the reciprocal of the input 
pulse duration. How big the ratio of these two quantities must be de- 
pends upon how much overlap there is in the output pulses and there- 
fore upon how much equalization is required to separate the pulses at 

baseband. 

We can model the complex amplitude of a given spatial mode at the 

guide input as follows: 

e iB (t) = V^KT)c(0. (1) 

In the above example, («?(£))* represents the modulation and c(t) is 
a complex Gaussian random process which represents the incoherent 
carrier. By definition of an incoherent carrier, we have 

(c(t)c(t + u)) = 0, (c(t)c*(t + w)) = R c (u). (2) 

The Fourier transform of R c (u) is what is called the incoherent source 
spectrum, shifted to baseband. 

The input complex amplitude of (1) produces an output from 
the guide in the corresponding mode having the following complex 
amplitude 

Cout(0 = J*in(t')h g (t - t')dt'. (3) 

In eq. (3), h g (u) is the guide bandpass impulse response for the mode 
under consideration. Equation (3) follows from the fact that the guide 
is linear in voltage. 

The average power at the guide input in the given mode is (averaging 
over the fluctuations in the incoherent carrier) 

(Pin(t)) ^ <ein(*)4(*)> = m(t)R e (0) = R c (0)lE a k h pln (t - kTn (4) 

In (4), the modulation m(t) is a sequence of nonoverlapping pulses 
(modulated on or off) where h Pin (t) is the pulse shape, a* assumes the 
value zero or one for each k, and T is the pulse spacing. The average 
power at the guide output is 

(Poutit)) = (jjtUt')h g (t - 04(O*!(* - t")dt'dt'\ ■ (5) 

It is reasonable to assume that the following approximation holds: 



<em(04(n> = {^W)<n<Mnc*{t")) = <W)<^n 

R c {t' - t") ttm{t')R c {t' - t"). (6) 
This approximation is valid since the coherence time of a typical 
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incoherent source such as a GaAs LED is of the order of 10 -13 seconds, 
while modulation pulse widths of interest here exceed 10 -9 seconds. 
Substituting (6) into (5) we obtain 

<7>out(0) = L a k h Po Jt - kT), (7) 

where 

h Poat (t) =fh Pin (t') \JRc(t' - t")h g (t - t')K(t - t")dt"~\ dt'. 

Thus, the average output power is a linearly filtered version of the 
average input power. That is, the average output power consists of a 
sequence of pulses which add linearly even if they overlap. 

Thus far, we have considered the average output power, averaging 
over the fluctuations in the incoherent source power output. We next 
consider the effect of those fluctuations on the equalized-detected 
current. We shall show that these source fluctuations will produce 
negligible deviations in the equalized-detected current from its mean 
provided that the source bandwidth is sufficiently large. 

We can write the power at the guide output in the mode under 
consideration as the sum of the average power of eq. (7) and a deviation 
from this average b(i): 

Po«t(0 = <Pout(0) + KO- (8) 

If this power falls upon a detector, it produces a current which is 
proportional to p ou t(0 (neglecting shot noise). This current will pass 
through a filter which performs an equalization function and/or band- 
limiting. The baseband filter output voltage will therefore be 

Vont(t) = zjp 0Xlt (t')h b (t - t')dt', (9) 

where hb(t — t') is the baseband filter impulse response and z is an 
arbitrary proportionality constant. 

The mean baseband output voltage is given by 

<*Vt(0> = zf(Pout(t'))h b (t - t')dt'. (10) 

The mean square deviation of the baseband voltage from its mean 
is given by 

fc(0> - [^out(O)] 2 = *\{t) ^ Z* f(Pout(t')Pont(t")) 

X h b (t - t')h b (t - t")dt'dt" - [>out(0>] 2 - (11) 

Thus, to calculate the ratio of the mean voltage to the rms deviation 
in order to determine whether or not the deviations are negligible, we 
need the correlation function of the power p ou t(0- 
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In order to calculate this correlation function, we must recall that 
c(t) denned in (1) is a complex Gaussian random process and satisfies 
[in addition to (2)] the following 

(c(t)c*(t')c(t")c*(t'")) 

= R c (t - t')R c (t" - t'") + R c (t - t'")R c (t" - t'). (12a) 

Using (1), (3), (4), and (12a) we obtain 

<Po«t(0Po«t(<')> = <Po»t(0)<Po»t(O> 



+ 



fm(a)R e (a - (3)h g (t - a)h*(t' - 0)dad(3 



(12b) 



To obtain some numerical results, we assume that the input pulses 
h pln (t) defined in (4) are Gaussian in shape and that the guide mode 
impulse response corresponds to that of a dispersive medium having a 
group delay to at the optical source center frequency and a dispersion 
7 2 within the optical band of the source. Further, we shall assume that 
the source spectrum is Gaussian in shape. That is, we shall assume the 
following : 

h PiB (t) = exp (-t 2 /2<r 2 ), 

(2.36<r = input pulse width between £ points) (13a) 

S[h.(i)} = H.( a ) 

= exp ( - j ^f ) exp ( - jr «) exp ( - ^B~ 2 /2) (13b) 

h '® = \o /-l-uw-a 6XP [ " (t ~ r "> 2 /20V + Br*)l n . , 

V2ir(^7 2 + B„ 2 (13c) 

R c (u) - exp (-w 2 £ 2 /2), 

where 

£F{ } = Fourier transform. 

In eq. (13), we have already assumed that the source bandwidth B s 
is much greater than the reciprocal of the input pulse width a. We also 
shall assume that the guide bandwidth B g is much larger than the 
source bandwidth B„. Since t„ represents an absolute propagation 
delay from input to output, we shall neglect it as irrelevant to the 
problem at hand. Therefore, the only significant parameter in the guide 
impulse response h g (t) is the dispersion y 2 . 

If we insert the particular functions of (13) into (7) and (12b) using 
(4), we obtain 

<Pout(0) = ^2 E a k exp |-i(«- kT) ViyJE + <7*]} (14a) 
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(i.e., a sequence of Gaussian-shaped on-off modulated pulses having 
width ^y*B 2 + o- 2 ), where z 2 = an arbitrary proportionality constant 
which we shall henceforth set to unity. 



<Po«t(0Pout(O> = [<Po«t(0>? 



1 + exp 



(t - t'YB 




(14b) 



Thus, h PoM (t) of (8) is in this case the Gaussian-shaped pulse in (14a). 
Looking at eqs. (13a) and (14a), we see that y*B 2 /a 2 « 1 implies that 
little pulse broadening has occurred in propagation; y*B 2 a /a 2 ~2> 1 
implies that considerable pulse broadening has occurred in propagation. 
Now using (14) in (11), we obtain 

— Z27a 7 <w at , (15) 

,w J Dw(03*i(< - 0*i(* - r)e-<''-'" J, « ,/a 

where 

u 2 = 2B 2 [y i B 2 /<r 2 + I]" 1 . 

Recall that (v ou t(0) is the average baseband (detected and equalized) 
voltage produced by the power output in the mode under consideration, 
and that it consists of a sum of on-off modulated pulses given in (14a). 
In addition, a 2 v (t) is the mean squared deviation of this baseband 
voltage from its mean. Thus, (15) is effectively a signal-to-noise ratio. 
For a typical broadband source and an equalizer having a bandwidth 
comparable to 1/ (pulse spacing-T), we can treat the Gaussian term of 
the integrand in the denominator of (15) as a delta function having 
area V27r/w. Then (15) becomes 



Mff . [j(pUi'))h(t-w]' 



w -&'■ (16) 



u 

We can evaluate (16) for particular equalizers, hb(t), and particular 
output power pulse widths (y*B 2 + a 2 ). We recognize from (16) in 
general that the equivalent noise b(t) of (8) which must be added to 
the average output power is a signal-dependent noise with correlation 
function 

(b(t)b(t')) 2 ^E t<p mt (t))yt(t - t>), (17) 
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where 

u 2 = 2B\I\_^B 2 J« 2 + 1]. 

If the power output pulses h PoM (t) overlap significantly, then this 
signal-dependent noise will become stationary when all the pulses are 
on, which should simplify the calculation of (16). At the other extreme, 
if the output pulses h Po „,(t) do not overlap (i.e., y A B 2 + a 2 « T 2 ) and if 
the baseband equalizer hb(t) is taken to be a matched filter (matched 
in shape to the output pulses), then the signal-to-noise ratio is given by 

C^#l 2 = V2„B. (18) 

[if T 2 » y*B 2 + a 2 , h b (t) = h PoM (t)J 

In general, for a given desired equalized pulse shape, the equivalent 
noise will be negligible if the product of the optical source bandwidth 
B. and the input pulse width a is sufficiently large. For practical cases 
of interest, this product is on the order of 10 4 to 10 8 . For reasonable 
amounts of equalization consistent with other noise considerations 
(shot noise, thermal noise, etc.), we can treat the guide as being linear 
in power even if the pulses overlap, i.e., we can neglect the equivalent 
noise b(t). 

When more than one mode is present, we simply add the individual 
mode output powers, since the detector current is proportional to the 
sum of the powers in all the modes. If the optical source is spatially 
incoherent, the equivalent noises b(t) in each mode may be uncorre- 
cted, t In that case, the requirements upon the product of the source 
bandwidth and the input pulse width are less stringent. This is par- 
ticularly true if the pulse spreading resulting from dispersion dominates 
the spreading resulting from the differences in group delay among the 
various modes. 

A simple interpretation which may prove useful follows. It was 
easy to obtain (7), which showed that the guide was linear in power 
if we averaged out the source fluctuations. Since the source is very 
broadband, we can think of it as a sum of independently fluctuating 
sources separated by a frequency spacing equal to the bandwidth of 
the modulation pulses at the input. Thus, the output power is the sum 
of the fluctuating powers associated with each equivalent indepen- 
dently fluctuating optical source. The average output powers associated 
with these equivalent sources add systematically, while the indepen- 

t That is, if the optical source is close to the fiber, each fiber mode effectively sees 
an independent carrier, and therefore we obtain averaging of the fluctuations in 
these carriers when the mode powers at the output of the guide are added. 
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dent fluctuations about the average add at random. Thus, for a 
sufficiently large number of equivalent sources — corresponding to a 
large product of optical source bandwidth and input pulse width — the 
total fluctuations become small compared to the average power. In 
effect, one has a frequency diversity system. 

IV. CASE III : MODE-LOCKED COHERENT SOURCE, MULTIMODE GUIDE, 
MODE COUPLING 

In this example, we consider a spatially coherent mode-locked laser 
source and a multimode guide. Unlike example I, we assume con- 
siderable mode coupling. Once again, the rationale of using a multimode 
guide with a single mode source may be the inability to stably match 
the source to a single-mode guide. Before proceeding, we must model 
the transmission properties of a multimode guide with random coupling. 

Very little is known about the complete statistical properties of the 
guide under consideration here. Any particular guide, which is linear in 
voltage (field), can be characterized as having a set of modes associated 
with an ideal guide having no geometry perturbations (which are the 
source of coupling). The input and output complex envelopes of the 
corresponding input and output optical fields can be expanded using 
the orthogonal guided modes, which together with the continuum of 
radiating modes form a complete orthonormal series. We can relate the 
complex amplitudes in each input and output mode by a matrix im- 
pulse response. That is, calling the complex amplitude in mode k at 
the input e* to (0, and the complex amplitude in mode j at the output 
«y.«i(0i we have 

e/„«(0 = Z j *k ia (t')h jk {t - t')dt', (19) 

where hju{t — t') is the bandpass impulse response from input mode k 
to output mode j. 

In order to proceed in the analysis to follow, we need at least the 
fourth-order joint statistics of the random processes h jk (t). (The reader 
is cautioned that the term random process refers to the fact that the 
actual hj k (t) for each j and k will be different for different guides be- 
cause of the random mode coupling. For any particular guide which 
does not change its physical parameters in time, &,*(<) is a fixed but 
a priori unpredictable function of time. All averaging and references to 
statistical properties refer to ensembles of guides whose gross physical 
properties are alike.) 

As mentioned above, little is known about the statistics of the 
hjk(t). Rowe and Young, 1 Personick, 2 and Marcuse 3 have shown in 
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various analyses that, for a particular j and k, the Fourier transform 
of hjk(t), Hj k (co) can be considered a stationary random process under 
various restricted conditions which include the assumption that the 
optical bandwidth being used is not too large. That is, one may argue 
under restricted conditions that 

(£*(«)£*(« + *)> = 0Ji(«O (20) 

(where the averaging is over an ensemble of guides having identical 
gross properties). 

In another analysis, Marcuse 4 has shown that, for a particular j and 
k and for a sufficiently long guide (so that enough mode coupling has 
taken place), one has 

(|ff ifc ( W )| 4 )w2(|H, fc («)|»>». 

This last result is consistent with (but certainly not a sufficient 
condition for) the possibility that H jk {u>) is a complex Gaussian random 
process. 

Based on this admittedly scanty evidence which should certainly be 
explored in more depth, we shall assume that the Fourier transforms 
Hjk(<a) of the hj k (t) satisfy the following conditions which would be 
satisfied if the Hjk(u>) were joint complex Gaussian random processes 

(HMHUa + a)) = 0, (H,- k (o>)Ifi m ( U + «r)> = S jk U<r) 

(?*(«)&!(« + <r)H no (o> + </)#;> + o> / / ( 21 ) 

= Sjklm((T)Sno P q(v" ~ O + Sjkpq((r")S n olm(<r ~ 0"). 

It is hoped that the results which we shall next derive will be qualita- 
tively valid for actual multimode guides with random mode coupling. 
The guide input optical field complex amplitude is given by 

6in(0 = e p (p)e T (t), (22) 

where e p (p) represents the spatial variation of the field over the guide 
input plane and e T (t) represents the time variation of the field and 
includes the modulation. We shall assume that the modulation con- 
sists of a sum of nonoverlapping on-off modulated pulses with spacing T 

e T (fl = L *Mt - kT), 
a,k = or 1 for each k. 

We can expand the input field in the guide modes as follows : 

ep(p) = E e*0*(p) + unguided remainder, (24) 

where </>*(/>) is guided mode k. 
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From (20) and (24), we obtain the complex amplitude in mode j at 
the guide output in terms of the input complex envelope and an im- 
pulse response 

*/.«(<) = £ [ t k e T (t')h jk (t - t')dt' 
k J 

- f«r(W - W, where h } {t) £ £ 6*M*)- (25) 

Since /iy(J) is a weighted sum of individual responses h jk (t), it follows 
from (21) and the linearity properties of the Fourier transform that the 
transform Hj(ia) of hj(t) must also satisfy 

(Hj(w)H]( u + <r)> = Sj(a), {HMHiit* + *)> - 
<ffi(«)ffj(« + *)#,(« + OA|(- + O) (26) 

= SMSfo" - a') + 8fa")Bfa ~ a'). 

We next show that the guide may be considered under restricted 
circumstances to be linear in power with a well-defined output pulse 
shape even though the output pulses overlap and even though there is 
unpredictable mode coupling. 

First, we can write down the power at the guide input and at the 
guide output in mode j: 

Pin(t) = \j\ c p(p) 1 2d l\ £ a * I M< -kT)\ s (27) 

(since hin(t)h* D (t - kT) = for k ^ 0, and a k = or 1). 

P*.t(0 = j it{t')tr{t")hi{t ~ t')h*(t ~ t")dt'dt". 

Next, we can ensemble average p/ o .,(0 over the ensemble of similar 
guides to find the average power response. 

In Reference 2, it is shown that (26) implies that 

(hj(t - t')h)(t ~ t")) = 8j (t - t')S(t' - t"), (28) 

where s,-(t) = &~ l {Sj(v)) (inverse Fourier transform). 
Using (28) and (27) we obtain 

(pa-(0> - £ «i**-« - kT), 

where 

hjjf) = f | h ia (t') \ 2 sj(t - t')dt'. (29) 

We thus see that the input power consists of a sum of on-off modu- 
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latcd overlapping pulses, and the average output power (averaged over 
an ensemble of guides with identical gross physical properties) consists 
of a similar sum of on-off modulated pulses which in general may 
overlap. Thus, in this ensemble-averaged sense, the guide is linear in 
power with impulse response Sj(t). We must next investigate how an 
individual guide in the ensemble can deviate from this average and 
under what conditions these deviations can be neglected for com- 
munications purposes. 

In order to study these deviations we must consider them in the 
context of a detector followed by an equalizing (or simply band- 
limiting) filter. Since the detector produces a current which is pro- 
portional to the linear sum of the powers in each guide mode at the 
output (neglecting shot noise), we consider the response to one mode 
only for the moment. The voltage at the equalizer output is related to 
the output power in mode j as follows: 

Vout(t) = zfp Jm At')hd*.-tmXt - t')dt', (30) 

where /tdet.-fiit.(0 is the detector-filter impulse response and z is an 
arbitrary constant. 

The average (over an ensemble of guides) voltage produced is a sum 
of on-off modulated pulses, since we have already shown that the 
average output power in mode j is a sum of on-off modulated pulses 

<tw(*)> = L a k h VoM (t - kT), (31) 

where 

huJ& = fh joat (nh det ..mt.(t - t')dt'. 

The mean squared deviation from this average voltage is given as 
follows : 

<*(«) - M»(4) - <*w(<)> 2 

= z 2 j j ' <put')pun)h A «..mtXt - n 

XAd*.-!iH.(< ~ t")dt'dt" - (tWt(0>*. (32) 

To calculate the mean squared deviation, we need the correlation 
function of the output power in mode j. From (27) we obtain 

(pUa)pU0)) = IJJjler(t')et(t")e T (t'")e:(t""n 

X Ma - t')h*j(a - t")hj(p - t'")h*AP ~ t""))dt'dt"dt'"d.t"". (33) 
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In order to obtain numerical results, we must make some assump- 
tions to facilitate the products and convolutions of (33). We assume 
that the guide power impulse response Sj(t) for mode j is Gaussian in 
shape. It has been shown 1 ' 2 that this should be the case under re- 
stricted conditions for long guides. We assume that the input pulses, 
hin(t), are Gaussian in shape with a width less than £ the pulse spacing, 
T, so that the previous assumption that they do not overlap is not 
violated for practical purposes. That is, we assume 



h in (t) = e-' 7 > 2 °\ 3<r < T = pulse spacing 
Sj {t) = g-*W 7 > 3a. 



(34) 



°- + y> 



(35) 



From (33), (26), and (34) we obtain 

<P*«(<0> = zZa k exp-Ua- kT) 2 Al" 
<Pio„.(a)P/o».(£)) = U 2 EZE CLka k - m aiai +m 

L kml 

f (a - kT) 2 } f (8 - IT) 2 

XeX P-{ 2(, 2 /2 + ^) j eXp -{ 2(;V2 + 7 2 ) . 

Xexp j 2(.V2 + tV 2 IJ 

+ (psMXpuM)- 

The approximations of (35) become equalities when the width of the 
average power impulse response y becomes large compared to the input 
pulse width, a, i.e., when there is a lot of pulse spreading in propagation. 
We shall soon see that this will be the case of interest in this example. 
Before attempting to use (35) in (32) to evaluate the magnitude of the 
deviations of a particular guides power from the ensemble average, we 
can make some simplifications and comments upon (35). Using the 
assumptions yy><r, we have 

<?/.»,(«)) = * Z a k exp { - (a - kT) 2 /2y 2 } (36) 

R><S*,0) £ k,(a)p,„(3)) - (PioMXpUP)) 

= z 2 ELL a*a*_ m a i a i+m et-c-* 7 ') 2 / 2 > , )ei-(/3-/r) 5 /2^j 

kml 

\/ g- (a-0-m D 2 / 2<r J 

If we assume that a is small compared to the pulse spacing T and if 
we assume that the detector-equalizer combination passes frequencies 
only up to the inverse of the pulse spacing T, then we can make the 
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further approximation 

Rb(a, 0) = z 2 EEL a k a k - m aiai +m 

klm 

X [ 6 ( a -tr)'/27V(Mr)'/27 J ] g ( Q _ _ mT)^2T<r 2 . (37) 

What we have shown so far is that the power in mode j at the guide 
output can be considered to be of the form 

pUt) = ((pUD) + Ht), (38) 

where 

(b(a)b(f})) = R b (a ) 0), 

where b(t) represents the deviations of the power in a particular guide 
in mode j from the ensemble average. 

Combining (36) and (38) with (32), we obtain the ratio of the 
(average voltage) 2 at the detector filter output to the mean square 
deviation from this average voltage 

8 a (tfout(O) 8 
N o*(t) 

\ /E a*e- ( ''- fcT)J/272 W-fut.(* - t')dt'Y 
" (V2Ta)/EEE a kai a k - m a l+m e^'-^'^e^"-^'^ 

X 8(t' - t" - mT)h det ..intXt ~ t')h det ..mt.(t - t")dt'dt". (39) 

It is clear that, whatever the equalizing filter is, this ratio increases 
with decreasing a. Thus, as the input power pulses to the guide become 
narrow, two things happen: the average output pulse widths become 
independent of the input pulse width and the deviations of the output 
power in mode j from the ensemble average become negligible. Exactly 
how small a has to be depends upon how much the average output 
power pulses overlap and how much equalization we are therefore 
using. In the extreme case of no output pulse overlap, assuming that 
the equalizer response /idet.-nit.tf) is matched to the average pulse 
power, we obtain 

T > 3t (no pulse overlap at output) 

4 = 3—- for a « 7 (output pulse width » input pulse width) (40) 
N a 

W-nit.(0 = ***(-*) (matched filter equalizer). 

Obviously, the guide acts linearly in power if the output pulses 
don't overlap, but (40) shows that the output pulses take on a well- 
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defined shape (i.e., the deviations from the ensemble average are 
negligible) in spite of the random mode coupling. For cases where the 
output pulses overlap, the conditions for the deviations from the 
ensemble average to be negligible are more stringent, i.e., the signal-to- 
noise ratio (39) is less than the special case (40). 

Summarizing, we started with a mode-locked laser putting out 
pulses which were much narrower than the spacing between them and 
on-off modulated. This optical field excited a guide with random mode 
coupling. We modeled the transfer function relating the input field 
complex envelope to the complex amplitude of a particular mode at the 
output as having specific properties (26) associated with a complex 
Gaussian random process. This model was justified only in the sense 
that it was consistent with available but scanty analytical results on 
guides with random coupling. We showed that the ensemble average 
output power in the mode under consideration looked like a linearly 
filtered version of the input power. That is, the average output power 
was a linear sum of pulses which could in general overlap. Thus, on the 
average, the guide looked linear in power for digital communication 
applications. We showed that the deviations in a particular guide from 
this ensemble average linearity behavior would be negligible provided 
the input pulses were very narrow compared to the width of the guide 
average power impulse response. How narrow the input pulses had to 
be depended upon how much equalization was required to separate the 
output pulses. t 

It is clear that, since the detector adds the powers in each output 
mode, the total power will be a linear sum of pulses if the individual 
mode powers are. In addition, we may suspect that the deviations from 
the ensemble average in each mode may add randomly while the 
average powers add systematically. Thus, some improvement in the 
signal to "noise" ratio may accrue from this spatial diversity. 

An interpretation of what is happening to make the deviations 
negligible is the following: Since the guide average power impulse 
response for output mode j, Sj(t), is the Fourier transform of the two- 
frequency correlation function Sj(a) defined in (26), we can interpret 
the reciprocal of the width of s 3 (t) as the bandwidth difference over 
which the guide transfer function between the input field and the 
output mode j becomes uncorrelated. When we use narrow input 
pulses compared to the width of the average power impulse response, 
we use a lot of bandwidth compared to this correlation bandwidth and 

> Remember that the output pulse shape becomes independent of the input pulse 
shape as the input pulses get narrow. 
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thus obtain frequency diversity. As the input pulses become very 
narrow, the output pulses become fixed in shape equal to s } (t), but the 
diversity keeps increasing, resulting in averaging out of the deviations 
from one guide to the next. 

V. CASE IV : INCOHERENT SOURCE, MULTIMODE FIBER WITH MODE 
COUPLING 

In this example, we consider an incoherent intensity-modulated 
source exciting a multimode fiber with random mode coupling. We 
assume that material dispersion is negligible. (Since we shall be con- 
sidering a wideband source, we may question the physical reality of 
neglecting material dispersion. On the other hand, the qualitative 
results to follow may provide insight into more general cases.) To 
simplify what will prove to be a somewhat complicated analysis, we 
shall consider the response in a particular guide mode, j, at the output, 
to the field in a particular mode, k, at the guide input. Extension to 
consideration of the total input field and the total response should be 
straightforward, using the techniques outlined below. 

The input field complex amplitude in mode k is of the form 



e ki M = <rn(fic{t), (41) 

where m(t) is the modulation and c(t) is the optical incoherent carrier 
which satisfies 

(c(t)c*(u)) = R c (t - u) (42a) 

(c(t)c*(u)c(t')c*(u')) 

= R c (t - u)R e (t' - u') + R c (t - u')R c (t' - u). (42b) 

The complex amplitude in mode j at the output due to the input 
field in mode k is given by 

«/•*(*) = / '«* ta (OM* - f)*', (43) 

where the impulse response coupling mode k to j is assumed to satisfy 
the same statistics as were outlined in the first paragraphs of Section 
IV, i.e., 

(i) H jk (u) = 5{h, k (t)} 

(it) (H jk {u)H* jt {u + c)) = S jk (a) 

(Hi) (H jk (o>)H%(u> + a)H Jk (a + *')#£(« + O) (44) 

= S jk (a)S ik (a" - a') + 8 jk {a") S jk {a - <r') 

(w) (h jk (t)h%(t')) = s jk (t)8(t - t'), 
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where 

S jk (t) = S-MS;*(a>)). 

The average input power (averaging over the source fluctuations) is 
(ffihjfi) = m(t)R e (Q) = [Z a h h iQ (t - kT)]R e (Q), (45) 

i.e., a sequence of on-off modulated pulses where a k = or 1, T = pulse 
spacing. 

The output power is given by 

PjoJt) = [Jj^J^(nc(t')c*(t")hj k (t - t')h* t (t - t")dt'dt". (46) 

If we average over the source fluctuations and the guide statistics 
we obtain, using (44) and (6) 

<PuM = Jm(t')R c (t' - t")(h jk (t - t')h* k (t - t"))dt"dt' 

= Jm(t')R c (0)s jk (t - t')dt' = £ a k h joM {t - kT), (47) 
where 

hj Jt) = R„(0) f hin(t')s 3k (t - t')dt'. 

Thus, on the average, the output power looks like a linearly filtered 
version of the input power with impulse response s> k (t). 

This average output power will produce an average detected- 
equalized voltage which is also a linearly filtered version of the input 
power. We next show that, if the optical bandwidth of the incoherent 
source is sufficiently large, then the deviations of the detected equalized 
voltage from its average, because of fluctuations in the source and 
deviations of the impulse response of a particular guide from the 
ensemble average, can be neglected. 

As before, the detected-equalized voltage is given by 

«W(0 = *Jp m t(t')h M ..mtXi ~ W, (48) 

where z is an arbitrary constant. 

The mean detected-equalized current and the mean square deviation 
from the mean are given by 

<»out(0) = £ a k h Voul (t - kT), 
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where 

h+Jfi = fhtJfihA^mtXt ~ t')di' 

<r 2 ,(t) = <«gu t (o> - <»out(o> 8 = f (punpun) 

X Adi..fm.(< - «0*d-.-fii».(* - t")dt'dt" - (v out (t)) 2 (49) 

To calculate the mean square deviation we need the correlation func- 
tion of pj wi (t). Using (6), (42b), and (44) we obtain the following com- 
plicated expression: 

(vut')vun) = bunxpun) [ i + |fic( ^(o) 01 ' ] 

+ /&"(*'»#(«' - »)•#(«" - r)dr j \R c (u)\Hu 

xf\Rc(u)\ 2 du. (50) 

In order to obtain numerical results we must make some assumptions 
as to the shapes of s ]k (t), h ia (t), and R e (t). We shall assume the 
following : 

fti»(0 = exp — t 2 /2a 2 (Gaussian-shaped input pulses) 
fffflcCOl = exp —o) 2 /2B 2 (Gaussian-shaped source spectrum) 
8/*(0 = exp — t 2 /2y 2 (Gaussian-shaped power impulse re- 
sponse — appropriate for long guides). 

Using (51) in (50) we obtain the expressions 

Ver 2 + 7 2 

(pUf)pun) 

= (pUt'))(pUt")){l + exp [-(«" - O 8 * 8 .]} 

-| Btay T T a k a j {e^ t "- t ' )2li - f2] e-^ k - j)T ^ lia> 

2\V 2 + t 2 

X e[~ 1/ ' a+ i ,2 ]i { ' + [('"-'')/ 2 ]-(*+^< r / 2 >l 2 + e-[(*-» r +(«"-'')] 2 / 4rf 

X e [-i/«r«+T*](«'+ 1 («"-«' )/2]-(fc+y)(r/2))2j > (52) 

Equation (52) is still fairly complicated, but we can make some 
general observations. The mean output power is proportional to the 
source bandwidth B„ and therefore the mean output power squared 
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will be proportional to B 2 . Looking at the mean squared output power, 
we see that the portions involving the double sum are proportional to 
B,. Those terms will become negligible compared to the mean output 
power squared as B. gets large. The term e xp —{t" — t') 2 B 2 s can be 
approximated by an impulse of area ^2-ir/B,. Thus, this term is also 
proportional to fi s _1 and will be negligible for large B,. Thus, for large 
B, we can approximate the mean square power by the mean power 
squared. Equivalently, a\{t) will become negligible compared to (v out (t)) 
as the source bandwidth becomes very large. Just how large the source 
bandwidth has to be depends upon how much overlap there is to 
average output pulses and therefore upon how much equalization has 
to be done. In the simple case where the output pulses do not overlap 
and where the equalizer is matched in shape to the output pulses, we 
obtain 

S (v out (t)) 2 V2B. 

V7 2 + e 2 o-V2 7 2 + «7 2 7 V2cr 2 + 7 2 
for 

T > 3Va 2 + 7 2 = 1.5 X output pulse width 
^dot.-fiit.(0 = hj mt (— t) (matched filter). 

From (53) we see that, in this special case, the ratio of the mean 
detected equalized (filtered) voltage squared to the mean squared 
deviation, because of source fluctuations and guide random coupling, 
will be greater than 1000 if the optical source bandwidth is more than 
1000 times the reciprocals of both the input pulse duration and the 
average power impulse response duration. 

Typically, the source bandwidth is more than 10" radians per 
seconds. Thus, the deviations are negligible for input and output pulse 
widths larger than 10 -10 seconds. Of course, more bandwidth is re- 
quired to make the deviations negligible when there is considerable 
pulse overlap. 

An interpretation is similar to previous cases. The requirement that 
the source bandwidth be large compared to the reciprocal of the 
duration of the input pulses allows averaging out of the fluctuations in 
the source. The requirement that the source bandwidth be large com- 
pared to the reciprocal of the guide average power impulse response 
duration gives the frequency diversity that averages out the deviations 
between guides resulting from random mode coupling. 
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VI. CONCLUSIONS 

We conclude that there are a number of interesting circumstances in 
which a fiber system, normally linear in voltage, can be considered 
linear in power for digital communication purposes. The input power 
consisting of a sum of on-off modulated pulses produces an output 
power (and thus a detected current) which is also a sequence of on-off 
modulated pulses, possibly overlapping. The output pulse shape is 
well-defined under the conditions described above, which usually 
amount to using enough optical bandwidth to have sufficient fre- 
quency diversity to average out source fluctuations and/or random 
mode coupling differences between guides. This power linearity allows 
baseband equalization (with the usual noise penalties) to allow use of 
the guide at higher bit rates than would be associated with the criterion 
that the output pulses must not overlap. Many assumptions made 
above may not be completely applicable to particular guides, but it is 
hoped that, qualitatively, some insight as to when power linearity may 
occur will be derived from these results. 
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